Transport properties of doped t-J ladders 
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Conductivity and Hall coefficient for various types of t-J ladders are calculated as a function of 
temperature and frequency by numerical diagonalization. A crossover from an incoherent to a 
coherent charge dynamics is found at a temperature T co h- There exists another crossover at Tpg 
below which a pseudogap opens in the optical spectra, induced by the opening of a spin gap. In 
the absence of the spin gap, T co h and the coherent weight are suppressed especially with increasing 
dimensionality. On the contrary, T co h is strongly enhanced by the pseudogap formation below Tpg, 
where the coherent Drude weight decreases with increasing dimensionality. The Hall coefficient 
shows a strong crossover at Tpg below which it has large amplitude for small doping concentration. 
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In high-T c cuprates, transport phenomena in the nor- 
mal state show various unusual properties in terms of 
standard normal metals. It is well known that the in- 
plane DC conductivity <t, is proportional to 1/T with T 
being the temperaturea, and the AC conductivity, has a 
broad incoherent part well fitted by <j(u>) oc 1/un. An- 
other peculiarity exists in the Hall coefficient, Rjj. At 
low T's, in spite of the electron density less than half 
filling, Rjj is usually positive, meaning hole-like carriers, 

fa large amplitude scaled as 1/6 with hole doping 
It also shows a crossover to a small amplitude at a 
rather high temperature, T cr , and T cr rapidly decreases 
with increasing Sou. These behaviors are incompatible 
with the standard Fermi liquid theory, implying that the 
carriers change their character with T due to strong cor- 
relations. The incoherent charge dynamics and strongly 
T-dependent Hall coefficient have been analyzed in terms 
of criticality near the Mott transition together with the 
effects of preformed pairsQ. Jaklic and PrelovsekH cal- 
culated <t(lu,T) for the two-dimensional (2D) t-J model 
using Lanczos diagonalization. Their results reproduced 
the above features of the high-T c cuprates, in contrast to 
the ID t-J model, which is known to have a very coherent 

Assaad and Imada showed a large crossover in the T- 
dependence of the high-frequency Hall coefficient, R* H , 
for the 2D iiubbard model using quantum Monte Carlo 
simulationslL9. At small 5's, R* H changes its sign twice 
with decreasing T from T=oo: from negative to positive 
at Tjj, and then back to negative at Taf- In the re- 
gion of Taf<T<Tjj, it is found that strong correlations 
substantially suppress a quantum coherence between dif- 
ferent spin configurations in the electron motion, leading 
to a large positive Hall coefficient expected for a spinless 
Fermion case. Within numerical results, the 2D Hub- 
bard model shows no definite indications for preformed 
pairing or superconductivity It neither reproduces the 
experimentally observed T-dependence of Rh- 

Recently several two-leg ladder compounds have been 
experimentally studied. A finite spin gap was con- 



firmed in the insulating phases, and the superconduc- 
tivity was reported for doped Sri4_ 2; Ca 2 ;Cu2404i under 
high pressureLU. The optical conductivity in the metallic 
region shows a similar incoherent character to the high- 
Tc 's while thejjesistivity follows p oc T 2 at low T's in the 
normal phasetj. 

In this paper, we study the conductivity and Hall co- 
efficient of the t-J two-leg ladders. The Hamiltonian of 
the t-J ladder under a magnetic field, B, perpendicular 
to the ladder plane is written as 



H = 



H.c. + J a Si ■ Sj 



(1) 



where the site sum is taken over nearest neighbor pairs, 
and the label, a=x(y), denotes the leg (rung) direction. 
The exchange constant is J x — J along legs and J y =J' 
along rungs. The hopping integral has a finite phase 
due to the magnetic field, Uj=t exp[— ie J } A(r) ■ dr], 
with the Landau gauge, A=B(— y, 0, 0). The periodic 
(open) boundary conditions are used along the leg (rung) 
direction. 

The conductivity at temperature T is defined by the 
Kubo formula as 



CT a p(uS) 



— {-{K a )5 a p 



-(n\J a \m)(m\J \n)), (2) 
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where V is the number of sites, Z— r Ire~ H ^ T , |n)'s are 
eigenstates of TL with energies E n , (K a ) is the thermal 
average of the kinetic energy along the a-direction, and 
J a is the paramagnetic current operator, rj is an adia- 
batic constant. The Hall coefficient is defined as 



Rh{uj) 



1 



B 0- xx (w)(Jyy{(jj) + <j2 {UJ) ' 



(3) 



Shastry et al. suggested its high-frequency lirpit as a mea- 
sure of T-dependcnt effective carrier densityLL3 
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We have numerically obtained all the eigenstates using 
the Householder method, and calculated the conductivity 
and the Hall coefficient using Eqs. The canonical 

ensemble is used to fix the electron number with all the 
subspaces of S* ot included. Spin gap, A s , is also calcu- 
lated as the difference of the ground state energy between 
the subspaces of S* ot — and 1. 

Let us start with the w-dependence of conductivity. 
Figure [j] shows <t xx (uj) at various T's for three typical 
cases: (a) uniformly coupled ladder (UL), J/t=J' /i=0.3, 
(b) rung-dominant ladder (RL) which has larger rung 
couplings than leg couplings, J 1 /t~2QJ/t—6.0, and (c) 
dimerized ladder (DL) where each chain is dimerized, 
J/t=J' /t=0.3, and J"/t=1.0. To characterize the dimcr- 
ization, we use the notation J=Jij<J" =Jij+i- At high 
T's above the temperature scale of A S) a xx (u>) has a 
broad incoherent part. At these temperatures, the con- 
ductivity is well fitted by the scaling function-proposed 
by Jaklic and Prelovsek for the 2D t-J modelli-3 



-C(w) 



(5) 



Here C(w) is a smooth universal function with a broad 
width, 7o of the order of several times of t. In 
our case, this may be well approximated by C(u>) ~ 
C (T)e- ( "/ 7 °) 2 [l- e -(To/T+T/ w )] for gmall v with a rather 

weak T dependence of Cq(T). The details will be dis- 
cussed elsewhere. The broad structure of C(u) is a con- 
sequence of rapid and incoherent charge dynamics, since 
as shown from Eq. (0), C(lo) is the current-current corre- 
lation function, C(J) = f dt{J(t)J(0))e iut . We note that 
this incoherent charge dynamics is consistent with the 
anomalous suppression of coherence near the Mott insula- 
tor derived in the scaling theoryO. It should also be noted 
that this scaling function with weak T, w dependence of 
C{ijj) corresponds to the relaxation rate proposed in the 
marginal Fermi liquid theoryO-J r _1 ~ max(w,T) at low 
to and T. 

Below a characteristic temperature, TpQ, the conduc- 
tivity <j xx (lo) turns to show a pseudogap at lo<Ll>pq as 
clearly shown in Fig. [I] (b) and (c). This is accompanied 
by the appearance of a prominent coherent peak in the 
low-o; region, particularly for Fig. [I] (b), and its weight 
grows with decreasing T. By changing J and J', we have 
found that Tp G ^A s ~w PG . Note that A s = 2.402 for 
(b) and A s = 0.812 for (c). The magnetic susceptibil- 
ity, x{T), is in fact substantially reduced at T < Tpq. 
In the case of Fig. [I] (a) , the pseudogap behavior is not 
prominent down to the lowest temperature used in the 
calculation because of the small spin gap A s /t « 0.002. 

Comparing Fig. [l] (b) (namely RL) and (c)(DL) at 
T < TpQ, it is noticed that the RL has a larger coherent 
weight. The observed l/(jj 2 -behavior of a xx (oj) at u/t<l 
in (b) is actually due to a finite r\. The ^-dependence 



at small rj indicates that this coherent peak is indeed 
a (5-function, ttD*S(uj), within our numerical results, al- 
though T > 0. In contrast, u xx (lo) for the DL in Fig. |l| 
(c) shows a larger incoherent weight, even below TpQ- 

The pseudogap behavior below Tpq indicates the for- 
mation of singlet bound pairs of electrons, which behave 
as singlet hard-core bosons. This is consistent with the 
exact diagonalization results at T = 0L3~ta. The part of 
lu<A s and T<Tpq is therefore described by an effective 
model of hard-core bosons with charge +2e, belonging to 
the universality class of Luther-Emery liquidslij. 

The bound pairs are formed on a rung in the RL's while 
on a dimerized strong bond in the DL's. Therefore it 
leads to different effective models, i.e., hard-core bosons 
on a single chain and on a ladder, respectively, with a 
small nearest-neighbor hopping t* and interactions V* as 
calculated in RefJlJ. The boson ladder and chain show an 
incoherent behavior, above a characteristic temperature, 
T co h, which is similar to the incoherent behavior of the 
Fermion ladder at T>Xpg discussed for Fig. |l| Whereas 
a dissipationless Drude peak, <r xx (uj)^irD* (T)5(ui) , ap- 
pears below T co h with Z?*(T)^const. T co h also corre- 
sponds to the crossover in the T-dependence of the ki- 
netic energy: roughly speaking, (K x ) ocl/T above T co h 
and (if^wconst. below it. Based on a finite-size scal- 
ing of (K X )(T), we found T co hC)Ci5 2 t*, consistent with the 
prediction of thg scaling theory for the metal-insulator 
transition in 1D0. 

The <t(u>, TLaf the boson chain was calculated by Zotos 
and PrelovsekEl They found that there remains a real 5- 
function of the Drude peak even at temperatures of order 
t* , and attributed it to the integrability of the model. 
Our results are consistent with their observation. 

We have calculated the conductivity by the same nu- 
merical method for hard-core boson chains and ladders. 
The results show that the Drude weight for the boson 
ladders is smaller than for the boson chains, suggesting 
a more diffusive character of the conductivity, whereas it 
turns out that T co h is higher for the ladders. The scaling 
theory for a single component systemQ predicted DotS 
and Tcohoc^ 2 in ID, and DccS and T co hoc<5 in 2D. The dif- 
ference between the results in Fig. [j] (b) and (c) implies 
a crossover from ID to 2D. 

Now let us discuss the Hall coefficient at w — > oo, 
R* H . Figure || shows the T-dependence of R* H for var- 
ious ladders including those used for Fig. [IJ. In the limit 
of T — > oo, R* H is finite and positive in all cases. For 
the t-J ladder with hole doping 5, it is indeed given by 
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R H 
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identical to the result for the 



n 

square latticetU, where it is positive if < 8 < |. With 
lowering T, R* H decreases and becomes negative at a cer- 
tain te™perature, Taf, as in the case of the 2D Hubbard 
modefcJ. It suggests the same origin of this behavior 
with the 2D case at T<Taf, i.e., the gradual formation 
of electron-like large "Fermi surface" in the presence of 
two spin species. 

In the UL's (J/t—J' /i=0.3 and 1.5), R* H remains neg- 
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ative at T<Taf- With decreasing T, it shows a minimum 
and finally converges to a large negative value at T=0. 
In the RL's (J/<=0.3, J'/t=6.0 and 20.0), R* H decreases 
with lowering T as in the previous case, but there appears 
a spike-like structure at around Tpq- This structure is 
more prominent for larger J"s. In the DL's, R H <0 in a 
very narrow T-region and R* H increases to positive and 
converges to a large positive value at T=0. 

The low-temperature behaviors, particularly at 
T<Tpg, depend on the coupling constants, meaning dis- 
tinct characters of charge carriers for the different cases. 
Since (— A" Q )>0 at any T's, the sign of R* H is determined 
by the sign of ([J x , J y ]), i.e., that of a xy . Considering the 
formation of bound hole pairs below Tpg, it is natural to 
assume that o~ xy in Eq. (|3j) is obtained as the sum of con- 
tributions from bound pairs and unpaired electrons. The 
T-dependence of R* H is therefore qualitatively explained 
from T-dependence of the density and mobility of these 
two types of carriers. 

As for the unpaired electrons, their density decreases 
exponentially at T<Tpq. Equation (^) shows their Hall 
mobility is determined by the processes of moving an un- 
paired electron first along one direction and then perpen- 
dicular to it. It is important to notice that the effective 
hopping matrix elements of unpaired electrons acquire an 
extra Fermion negative sign when they undergo a recom- 
bination scattering with a bound pair. Therefore, the 
Hall mobility of unpaired electrons is positive at T>Taf 
and T<S^Tpq while negative at Tpq<T<Taf, changing 
its sign twice. 

When the bound pairs are concerned, we recall that 
a xy at T <C Tpg is mainly determined by this part, since 
only few unpaired electrons are thermally excited. The 
bound pairs are formed on the rungs in the RL's, whereas 
within the legs for the DL's. The opposite sign of the 
Hall mobility at T < Tpg between these two cases may 
be attributed to these different local configurations of 
bound pairs, which will be discussed elsewhere. 

The Hall coefficient for small ui is shown in Fig. || for 
the same ladders in Fig. ^. The T-dependence is qual- 
itatively similar to R* H for the UL's and the DL's. On 
the other hand, the RL's have a different T-dependence 
at small oj's. This may be due to finite size corrections, 
since Rjj at small w's is very sensitive to the value of 
the adiabatic constant r\. Absolute value of Rh in Fig. || 
(a)-(d) also changes considerably with varying r\. 

Figure |J shows the hole doping dependence of R* H 
at T=0 calculated by Lanczos diagonalization for the 
UL's. When the electron density is small, R* H < 
0, meaning electron-like. When d > |, the value 
of R* H is close the classical value, i.e., R H ass =— yz^< 
R H (T=0)<R* H (T=oo)<0, showing little J-dependence. 
The Hall coefficient in this region is mainly determined 
by the electron density. On the other hand, at small S's, 
R* H shows a substantial J-dependence, and when J/t>l 
it is rather close to 1/6 except for its negative sign. This 
means the importance of electron correlation effects and 



is consistent with the previous conclusion that the bound 
pairs are effective charge carriers. A comment is neces- 
sary for the positive sign of R* H at 5=i and J/t<0.2. 
This may be due— to strong ferromagnetic fluctuations. 
As shown in RefO, the Nagaoka ferromagnetic ground 
state appears near half filling at small J's. The "Fermi 
volume" of the electrons with the majority spin is dou- 
bled, leading to a hole- like R* H . 

We finally discuss on how the Hall coefficient is affected 
when ladders are coupled. Figure | shows R* H {T=0) for 
various coupled ladders calculated by Lanczos diagonal- 
ization using Eq. (|J). Two two-leg ladders are coupled 
by interladder hopping and exchange Jl, to form a 
4x4-site cluster. The finite coupling of the ladders shows 
shifts of R* H toward the positive direction in general. 
A typical case is weakly coupled DL's and its equiva- 
lent configurations (for example, J"> J= J'=10Jl=0.3£ 
or 0.03T> J=J"=10J L =Q.3t or 0.6 J'> J= J"= J L =0M). 
A few exceptions with the negative R* H are the cases of (i) 
weakly coupled RL's, (ii) the configuration of weakly cou- 
pled plaquettes (e.g., the coupled DL at Jl<J<J"<J')- 
The crossover from a small to a large positive Rh at Tpc 
in 2D configurations is reminescient of the crossover in 
the cuprates at T = T cr , suggesting that it share*, the 
same mechanism due to preformed pairs below Tpqu. 

In summary, the frequency-dependent conductivity 
and Hall coefficient of various types of doped ladders 
were calculated by numerical diagonalization. There is 
a crossover from high-T incoherent to low-T coherent 
dynamics with the appearance of the Drude peak. Its 
crossover temperature T co h is enhanced by the spin gap 
and resultant pseudogap formations. T co h increases with 
increasing dimensionality when the spin gap is formed, 
although the coherent weight itself decreases. The Hall 
coefficient shows a strong crossover to large amplitudes 
at low temperatures below the spin-gap temperature. 
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FIG. 4. Electron-density dependence of R* H for the t-J 
ladder with uniform couplings, J=J'. 2x8 sites. 



FIG. 1. Frequency dependence of the conductivity of the 
t-J ladder for various choices of parameters, (a) J—J'=0.3t, 
2 holes in 2x6 sites, (b) J=0.3i, J'=6.0i, 2 holes in 2x6 sites, 
and (c) J=J'=0.3i, J"=t, 2 holes in 2x4 sites. The adiabatic 
constant is (a)-(b):r?=0. It, (c): O.Olt. (d) J(u) = f£ C(w')dw' 
calculated for (a). 



FIG. 2. Temperature dependence of the high-frequency 
Hall coefficient of the t-J ladder. The system size is (a)-(b): 
2 holes in 2x6 sites, and (c) 2 holes in 2x4 sites. 



FIG. 5. Hall coefficient of coupled t-J ladders in the cj=oo 
limit, K* H . Two ladders with 4x2 sites are coupled by inter- 
ladder hopping tL and exchange Jl- 



FIG. 3. Temperature dependence of the low-frequency Hall 
coefficient of the t-J ladder. The system size is (a)-(b) : 2 holes 
in 2x6 sites (5 = |), and (e) 2 holes in 2x4 sites (S = \). 
The adiabatic constant is r\ = O.lt. 
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oo (a) J=J'=at 

□a (b) J=J'=0.3t, J"=at 

( C ) J=J'=10J L =at, t L =0. 
um (d) J=J'=0.3t, J"=at 
^/t^/J^.l 
( e ) J=J"=0.3t, J'=at 
^1=3^=0.1 
♦♦(f) j=j"=J L =0.3t 
j'=at, t L =0.1t 
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